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The purpose of this paper is twofold. Firstly, we study the dynamical density fluctua- 
tions around the critical point(CP) of Quantum Chromodynamics(QCD) using dissipative 
relativistic fluid dynamics in which the coupling of the density fluctuations to those of other 
conserved quantities is taken into account. We show that the sound mode which is directly 
coupled to the mechanical density fluctuation is attenuated and in turn the thermal mode 
becomes the genuine soft mode at the QCD CP. We give a speculation on the possible fate 
of a Mach cone in the vicinity of the QCD CP as a signal of the existence of the CP on the 
basis of the above findings. Secondly, we clarify that the so called first-order relativistic fluid 
dynamic equations have generically no problem to describe fluid dynamic phenomena with 
long wave lengths contrary to a naive suspect whereas even Israel-Stewart equation, a pop- 
ular second-order equation, may not describe the hydrodynamic mode in general depending 
TT^ on the value of the relaxation time. 
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The analysis 2 -* of the experimental data 3 ) obtained in Relativistic Heavy Ion 
Collider (RHIC) at Brookhaven National Laboratory(BNL) has suggested that the 
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matter created by the relativistic heavy ion collisions is well described as almost ideal 



relativistic fluids with the ratio of shear viscosity rj to entropy density s being tiny. 
This discovery has prompted great interest in the origin and the true value of the 
viscosities including the bulk viscosity £ and so on, and people are now enthusiastic 
in elucidating dissipative effects 4 -* and thus in constructing a transport theory for 
strongly coupled systems. 5 )~ 7 ) 

A unique feature of the phase diagram of Quantum Chromodynamics(QCD) 
is the existence of a critical point(CP), which is yet mainly based on some effective 
models 8 ^ and a few lattice studies; 9 ) see Refs.lO),ll) and 12) for possible variants and 
alternatives. At the QCD CP, the first order phase transition terminates and turns to 
a second order phase transition. Around a critical point of a second order transition, 
we can expect large fluctuations of various quantities, and more importantly there 
should exist a soft mode associated to the CP of second order. 13 ) What is, however, 
the soft mode of the QCD CP? It has been established now that the QCD CP 
belongs to the same universality class as the liquid-gas CP and shown that the 
density fluctuating mode and generically hydrodynamic modes coupled to conserved 
quantities in the space-like region are a softening mode at the CP: 14 )' 15 ) The would- 
be soft mode, the a mode, is coupled to the density fluctuation 16 ) and becomes a 
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slaving mode of the density variable; 14 )' 15 )' 17 ) see Ref. 18) for another argument on 
the fate of the a mode around the CP. Furthermore it is also suggested 19 )' 20 ) that 
the bulk viscosity should show a divergent behavior around the QCD CP.*) 

In this paper, we shall show that the ultimate soft mode at the QCD CP may 
not be a sound mode that is directly related to the dynamical density fluctuation, 
and that the possible divergent behavior of the viscosities might not be observed 
through the density fluctuations; the sound modes are attenuated around the CP 
and would eventually die out at the CP; in turn, the diffusive thermal mode that is 
coupled to the sound mode would become the soft mode at the QCD CP. We should 
mention that Fujii and Ohtani 14 ) emphasized that the thermal mode as well as the 
density oscillation play the role of the soft modes in the QCD CP on the basis of 
the time-dependent Ginzburg-Landau formalism. However, the relative importance 
of these hydrodynamic modes around the QCD CP was not elucidated. 

Our investigation is based on an explicit use of the relativistic fluid dynamic 
equations for a viscous fluid. We shall show that even the so called first-order 
relativistic fluid dynamic equations have generically no problem to describe fluid 
dynamical phenomena with long wave lengths contrary to a naive suspect. 

In the case of nonrelativistic fluids, the nature of the critical dynamics around 
the critical point is rather thoroughly elucidated. 22 )~ 24 ) The static density fluctuation 
induces a strong light scattering in the vicinity of the critical point, which is called 
critical opalescence. The dynamical density fluctuations are investigated on the 
basis of the (non-relativistic) Navier-Stokes equation. 22 )' 24 )~ 26 ) Note that the fluid 
dynamics is a dynamical thermodynamics, which tells us that the density fluctuation 
is in general coupled to the entropy fluctuations creating diffusive thermal mode. 

We shall explore how the singularities of the thermodynamic values as well as 
the transport coefficients affect the dynamical density fluctuations around the QCD 
CP using relativistic fluid dynamics for a viscous fluid. We examine the spectral 
function of the density fluctuation in the equilibrium. Our analysis is actually an 
extension of that made for non-relativistic case 25 )' 26 ) to the relativistic case. We 
are not aware of such an analysis so far except for a simple estimate 27 ) using Euler 
equation without dissipation in a astronomical context. The use of the relativistic 
fluid dynamics for describing the dynamics around QCD CP should be interesting in 
view of the success of the fluid dynamics for RHIC phenomenology and the nature of 
the QCD CP mentioned above. 28 ) Admittedly, the analysis based on fluid dynamics 
is only valid in the fluid dynamical regime: k£ « 1, where k and £ denote the 
typical wave number of the fluid and the correlation length; that is, our discussions 
on the critical behavior are based on extrapolation from the fluid dynamical to the 
critical region: k£ >> l. 29 ) This extrapolation is known to make good sense for the 
non-relativistic case. 22 ) 

A set of the dynamical variables is thus the density n, temperature T and the 
fluid velocity u 11 (u^u^ = 1; fi = 0, 1, , , 3), while the sigma mode is not included as 
a dynamical variable because it is a slaving mode of the density. 14 )' 15 ) We shall see 
that the coupling of the density fluctuation with that of the entropy gives rise to an 
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important effect and is found essential for the description of the dynamics around 
the CP: It will turn out that the relativistic effects on the spectral function of the 
density fluctuation appears only in the width of the Lorentzian peaks due to the 
sound and thermal modes through the modification of the transport coefficients due 
to relativistic effects. Moreover, our detailed analysis of the critical behavior of the 
transport coefficients as well as the thermodynamic quantities will show that the 
possible singular behavior of the transport coefficients with relativistic effects turns 
out to be masked by the more singular behavior of the specific heats; hence we find 
that the sound mode will be attenuated , and only the Rayleigh peak due to the 
thermal fluctuation stays out around the QCD CP, which is found, to our surprise, 
precisely the same in effect as in the nonrelativistic case. 22 ) 

The forms of the dissipative relativistic fluid dynamic equations are far from be- 
ing established, although a great development has been seen in a couple of years. 5 )' 7 ) 
Relativistic fluid dynamic equations are the balance equations for energy-momentum 
and particle number, 

a M T^ = o, (m) 

d^N* = 0, (1-2) 

where T^ v is the energy- momentum tensor and iV M the particle current. They are 
expressed as 

T» v = (e + P)u»u v - Pg» v + r^, (1-3) 
= nu^ + i/", (1-4) 

where e is the energy density, P the pressure, the flow velocity, and n the particle 
density, the dissipative part of the energy-momentum tensor and the particle current 
are denoted by r^ v and v^, respectively. 

The so called first order equations such as Landau 30 -* and Eckart 31 ) equations are 
parabolic and formally violates the causality, and hence called acausal. Moreover, 
the Eckart equation which is defined for the particle frame where the particle current 
does not have a dissipative part shows a pathological property that the fluctuations 
around the thermal equilibrium is unstable, 32 ) while the Landau equation defined 
for the energy frame does not show such a pathological behavior. The causality 
problem is circumvented in the Israel-Stewart equation, 33 ) which is a second-order 
equation with relaxation times incorporated. One should, however, note that the 
problem of the causality is only encountered when one tries to describe phenomena 
with fast velocities which is beyond the scope of the fluid dynamics; nobody would 
actually expect that the fluid dynamics should describe phenomena with a velocity 
comparable with or greater than the mean velocity v of the constituent particles, 
though v is still less than the light velocity. The fluid dynamic equations describe 
the behaviors of the conserved quantities in even larger scale than the mean free 
path. Thus the phenomena which the fluid dynamics should describe are slowly 
varying ones with the wave lengths much larger than the mean free path. The 
problem of causality always takes place in short wavelength region such as shock 
wave phenomena, 34 -* so formally acausal fluid dynamic equations suffice in describing 
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fluid dynamical phenomena. In fact we will see that the results for fluid dynamical 
modes with long wave lengths are qualitatively the same irrespective whether the 
second-order or first-order equations are used. 

As for the instability seen in the Eckart equation, a new first-order equation 
in the particle frame constructed by Tsumura, Kunihiro and Ohnishi (TKO) 5 ) has 
no such a pathological behavior. 35 ) We employ Landau, 30 ) Eckart 31 ) and Israel- 
Stewart (I- S) 33 ) equation as typical equations, and TKO equation. 

This paper is organized as follows. In Sec. 2, we calculate the spectral function of 
the dynamical density fluctuation around thermal equilibrium state. We show that 
relativistic effects on the spectral function of the density fluctuation appears only 
in the width of sound modes, irrespective of the choice of the frame. In Sec. 3 we 
analyze the behavior of the spectral function around the QCD CP with the scaling 
laws, and show that the sound mode is attenuated while the thermal mode which 
is coupled to the density fluctuation stands out in the critical region. Section 4 
is devoted to discussions in which the fundamental reason is given why the sound 
mode is attenuated in the critical region in terms of the correlation length £ which 
diverges at the critical point. We furthermore suggest that the attenuation of the 
sound mode can be used to identify the existence of the QCD CP by the relativistic 
heavy-ion collisions; a possible suppression or disappearance of a Mach cone can be 
such an example. The final section is for summary and concluding remarks. We 
give a detailed derivation of thermodynamic relations in Appendix A, which are 
used in Sec. 2. In Appendix B, we show explicitly how the pathological behavior in 
Eckart equation affects the procedure of the calculation of the spectral function of 
the density fluctuations for an instructive purpose. In Appendix C and D, we present 
the detailed derivation of the spectral function in TKO equation in the particle frame 
and Israel-Stewart(I-S) equation, respectively. 



In this section, we derive the spectral function of the dynamical density fluc- 
tuations using some typical relativistic fluid dynamic equations for a viscous fluid, 
and discuss relativistic effects on and the frame dependence of the spectral function. 
A detailed derivation of the spectral function is given for the Landau equation in 
the first subsection. Leaving the similar detailed derivations to Appendix C and D, 
we present the results for the Tsumura-Kunihiro-Ohnishi(TKO) equation 5 ) and the 
Israel-Stewart equation 33 ) in particle frame in the subsequent subsections. 

2.1. In the case of Landau equation (energy frame) 

In Landau equation, the dissipative terms are given by 



§2. Analysis of dynamical density fluctuation 




(2-1) 



(2-2) 
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where rj is the shear viscosity, ( the bulk viscosity, k the thermal conductivity, 
T the temperature, fj, the chemical potential, and w = e + P the enthalpy density. 
A^ v = g^ u — u^n" is the projection operator to the space-like vector and = A^ v d v 
the space- like derivative (gradient operator). 

We calculate the spectral function of the dynamical density fluctuation around 
the thermal equilibrium state by the linear approximation for the deviation from 
the equilibrium. The following calculational procedure is an extension of the non- 
relativistic case described in the text book Ref.26). 

Let us write n(x) = no+5n(x), e(x) = eo+5e(x), P(x) = Po+5P(x), fi(x) = fio+ 
5(j,(x), and u^^x) = Uq + 5u^(x), where the respective quantities in the equilibrium 
state are denoted with q. For simplicity, let the equilibrium state be the rest frame 
of the fluid, Uq = (1,0); then owing to the relation u^Su^ = 0, it is found that Su^ 
takes the form 5u^(x) = (0, Sv(x)) with Sv(x) to be determined together with other 
quantities like dn etc. Then Landau equation Eqs.(lT) and (1-2) are reduced to 

96n + n V • 6v = K ^[^(5P) - V 2 (ST)], (2-3) 



dt wo wo 

01 



- vV 2 5v - (C + ^r?)V(V • 5v) + V{5P) = 0, (2-4) 



-o^ + -V 2 (<5P)-fv 2 (<5T) = 0, (2-5) 
at wo To 

in the linear approximation. Here we have used the following thermodynamic rela- 
tions. 

de = Tod{ns) + fiodn, (2'6) 
dP = noSodT + nod/j,, (2'7) 
w = T n s + n fj, , (2-8) 

where so is the entropy per the particle number in the thermal equilibrium. Note 
that the right hand side of Eq. (2-3) and the second term in Eq. (2-5) represent 
relativistic effects which are absent for the non-relativistic Navier-Stokes equation. 
In addition, the coefficient of the first term in Eq. (2-4) wo is replaced by the mass 
density po for the Navier-Stokes equation. 

Now we have five equations for seven unknown quantities, 5n, 8T, SP, 5s, and 
5v. In order to solve these equations, the thermodynamic quantities 5n, 5T, 5P 
and 5s are interrelated using thermodynamics. Choosing 5n and 6T as independent 
variables, we have 

5P(x) = ^5n(x) + ^^5T(x), (2-9) 

2 

Ss(x) = - W0 °f P 5n{x) + —5T(x), (2-10) 

where the following thermodynamic identities are used, 

( &P\ _ w oc 2 s ( dP\ _ woc 2 s ctp 



\dn J T n 7 \dT J n 7 
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where c n = To(ds/dT) n and cp = To(ds/dT)p are the specific heats at constant den- 

1 /2 

sity and pressure, respectively, c s = (dP/de) s ' the sound velocity, ap = — (l/no)(dn/dT)p 
the thermal expansivity at constant pressure, and 7 = cp/c n the ratio of the specific 
heats. Then Eqs. (2-3)-(2-5) take the form 

_ K M V 2 }5n + noV .Sv + k^-(1- C ^^)VHT = 0, (2-12) 
ot wqj wq 7 

w ^ - V V 2 5v - ( C + kv(V • Sv) + ^V<5n + ^^WT = 0, (2-13) 
ot 3 no7 7 

1 n 7 & ™o7 1 T ^ V 7 T ; ; V ; 

To utilize the correlations in the initial time to obtain the correlation function 
at later time t, we perform Fourier-Laplace transformation 

dr / dte- zt - ik - r 5n(r,t),---etc. 
-00 Jo 

Then we find 

(z + k 2 n^)5n(k, z) + inok ■ 5v{k, z) + k 2 K—( ^^ - 1)ST 
wqj wo 7 

= 6n(k,t = 0), (2-15) 

zwoSv(k, z) + k 2 r/5v(k, z) + (( + -rj)k(k ■ 5v(k, z)) + ik W °° s Sn(k, z) 

3 n 7 
2 

+ik W0 ° sap 5T(k, z) = w 5v(k, t = 0), (2-16) 

7 

wo^op _ k , K ± )5n{k ^ z) + { ^ _ fe2 c|op _ 1 
n 7 n 7 T 7 T 

= -^I^«5n(fe, t = 0) + ^5T(fc, t = 0). (2-17) 
n 7 T 

It is convenient to divide the velocity into longitudinal and transverse compo- 
nents 

5v(k, z) = bv\\ (k, z)k + Sv±(k, z). (2-18) 
The transverse component of Eqs.(2-15)-(2-17) reads 

zwo8v±(k, z) + k 2 r]5v±(k, z) = wo5v±(k, t = 0). (2-19) 

We shall not treat this equation, which admits a diffusive solution, in the present 
work since this equation is decoupled to the density fluctuation. We note that a 
complete treatment based on the mode-mode coupling theory 22 ) -24 ) of the critical 
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dynamics in the close vicinity of the CP involves the diffusive transverse mode as 
well as the thermal one which is also diffusive as we will see. 

The longitudinal component of Eqs.(2-15)-(2-17) can be cast into the following 
matrix form 



(5n(k, z) 
5v\\(k, z) 
ST(k,z) 

where the matrix A is defined by 



"07 



5n(k,0) \ 
8v\\(k,0) 



(2-20) 



A = 



( Z + k 2 K^ 

ik-^ 3 - 

"07 



. npc n I 
\w T I 



ikriQ 
z + vik 2 




apcjTo 



wo v 7 
ik apC " 



) \ 



7 



npc n 
uiqTq 



[z + /c 2 7x(l - 
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V 

(2-21) 

Here, we have introduced the longitudinal kinetic viscosity ui, and the thermal dif- 
fusivity x, 

A K 

(2-22) 



n = (C + ^ri)/w , 



X 



n c P 



Multiplying the inverse A 1 of the matrix A from the left in Eq.(2-20), we obtain 
the Fourier-Laplace coefficient of the density fluctuation 

2 

6n(k,z) = {(A-^n - ^( j4 - 1 ) 13 }^(fc, 0) + (A' 1 )^ (k, 0) 

+ ^(A- 1 ) 13 5T(k,0). (2-23) 



T w 

Now, we are interested in the spectral function of the dynamical density fluctu- 
ation 

S nn (k, u) = (5n(k, u)Sn(k, t = 0)), (2-24) 
where Sn(k,uj) is the Fourier transformation of the density fluctuation defined by 

dr / dte- iujt - ik - r 5n{r,t), (2-25) 

-OO J —CO 

and ( ) denotes the thermal average in the equilibrium. We see that S nn (k,uj) is 
obtained by taking the thermal average of Eq.(2-23) multiplied by Sn(k, 0). However, 
note that ST and 5n are statistically independent in fluids system which is established 
in Einstein fluctuation theory, 26 ) and accordingly 



Similarly, 



(ST(k,0)Sn(k,0)) = 0, 
(6vn(k,0)6n(k,0)) = 0. 



(2-26) 
(2-27) 
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Hence, we have 

S nn (k,v) = {(A~ 1 )n - ^(A' 1 ) 13 }{dn(k,0)5n(k,0)). (2-28) 

n 7 

The needed matrix elements (^4 _1 )n and (A^ 1 )^ may be calculated using the 
simple formula (A~ 1 ) n = 3^(^22^33 - ^23^32) and (A _1 )i 3 = 3^0(^12^23 - 
^-13^422)- Here, det^l reads 

det A = ^[z 3 + z 2 k 2 { 1X + n + - 2 X cWT } 

w T w 

+ zk 2 c 2 + k 4 c 2 X + ---], (2-29) 

where we have used the thermodynamic identity (A- 14), and '• • • ' denotes the higher 
order terms in k which we assume small because we are interested in the fluid dy- 
namical modes, det A can be nicely factorized in this approximation, 

det a ~ ^^(z + r R k 2 )(z + r B k 2 + ic s k)(z + r B k 2 - ic s k), (2-30) 

where 

r R = X, (2-31) 

= \\xil - 1) + n + c 2 s T (k/w q - 2 X a P )]. (2-32) 

Then we can write the Fourier-Laplace coefficient of the density fluctuation to 
second order in k 

5n(k,z) (z + r B k 2 + ic s k)(z + r B k 2 - ic s k) +zk 2 ^- fc 2 | 
5n(k,0) (z + r R k 2 )(z + r B k 2 + ic s k)(z + r B k 2 - ic s k) 

Performing the inverse Laplace transformation 

8n(k,t)/6n(k,0) = / dze zt 5n(k, z)/6n(k, 0), (2-33) 

2vTi Js-ioo 

we obtain the dynamical density fluctuation at t 

Sn{k, t) /5n(k, 0) ~ (1 - -)e- r * k ' H + - cos{c s kt)e- r * kH . (2-34) 

7 7 

Here, we have retained only the terms in the amplitudes to zeroth order in k. Since 
Eq.(2-34) is the density fluctuation in a stationary process, we can replace the time 
t by \t\. Therefore the Fourier transformation of Eq.(2-34) is given by 

Sn(k,u)/Sn(k,0) = (1 - 1 ^ 



7 7 cj 2 + r 2 k 4 



h | 1 (2-35) 

7 l (uj-c s k) 2 + r 2 k 4 (oo + c s k) 2 + r 2 ^ 1 ' [ } 
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Thus we finally obtain the spectral function of the density fluctuation 

S nn (k,u) = {5n(k,uj)5n(k,t = 0)} 

= <(<MM = o)) 2 >[(i-^ 21 Hk ' 2 



1 { FBk2 FBk2 
+ ^ { (oj-c s k)2 + r£k 4 + ( w + c s ky + JfF 1 v { ' 

We see that the spectral function have three peaks at frequencies uj = and uj = 
±c s k: The peak at uo = corresponds to thermally induced density fluctuations. This 
peak is called the Rayleigh peak, while the two side peaks at uj = ±c s k correspond 
to mechanically induced density fluctuation, i.e. sound waves. These two peaks are 
called Brillioun peaks. 

Now let us compare this results with that in the nonrelativistic case; 25 )' 26 ) 



S^(k,uj) = ((5n(k,t = 0)y)[(l--) 



7 uj 2 + x 2 k 4 

i r NRp r NRr,2 

+ -{ ^ + 5— } , 2-37 

7 (uj - c s kf + /^V (uj + c s fc) 2 + k 4. 



where 



^B NR =^[X(7-1) + ^ NR ], (2-38) 
^ NR =(C + ^)/p . (2-39) 
We see that relativistic effects appear only in the width of the Brillouin peaks: 

r B = r£ lR + jr£ a = r£ a , (2-40) 

where 

^ R ^[X(7-1) + ^], (2-41) 

with the longitudinal kinetic viscosity v\ defined in (2-22), and 

^La = I c 2 t ( k/u , _ 2 X a P ). (2-42) 

Firstly, the longitudinal kinetic viscosity is expressed in terms of the enthalpy density 
-wo in the relativistic case in place of the mass density p$. We call this modification 
the minimal relativistic (MR) effect. The other is a genuine relativistic effect 
which is absent in the non-relativistic case. This part comes from the right hand side 
of Eq. (2-3) and the second term in Eq. (2-5) which represent relativistic effects and 
vanishes if we take the light speed c — > oo. We see that the width of the Rayleigh 
peak is the same as the non-relativistic case 25 -*' 26 ) and the relativistic effects on the 
width of the sound modes tend to cancel with each other, so even the relativistic 
effect on the Brillouin peaks may be moderate. 
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Fig. 1. The spectral function in Landau and the minimal relativistic equation are shown by the 
solid line and the dashed line, respectively. The parameters are k = 0.1[l/fm], /tto = 200 [Me V], 
To = 200 [Me V], f?/(noSo) = C/( n <) s o) = 0.3 and «To/(noSo) = 0.6. Relativistic effects does not 
show up in the Rayleigh peak but enhance the Brillouin peaks. 



In order to see the relativistic effects 5T]f - quantitatively, we calculate thermo- 
dynamical values, c s ,ap and 7 by the equation of state (EoS) of massless classical 
ideal gas, e = 3P = 3nT, so we have c p = 4, c n = 3, ap = 1/Tq, c s = y/l/3 and 
so = 4 — hq/Tq . Then the widths in the minimal relativistic case are given by 

K 

4n ' 
1 



nMR _ 
1 B — 



(4r/ + 3C + KT ), 



(2-43) 
(2-44) 



(2-45) 



24n T 

and the relativistic corrections are found to be 

Srt = \c 2 s T ( K /w - 2 X a P ) = 

We see that the relativistic effects in the Brillouin peaks is solely due to the thermal 
conductivity and has a negative value, implying that the Brillouin peaks get enhanced 
with a smaller width by the relativistic effect. It addition, this relativistic term 
exactly cancels out with the minimal relativistic part due to the thermal conductivity 
in Eq.(2-44); thus we see that the thermal conductivity does not affect the Brillouin 
peaks in net for the massless classical ideal gas. 

Figure 1 shows the spectral function Eq.(2-36) and the minimal relativistic case 
with the parameter set k = 0.1[l/fm], ^0 = 200[MeV], T = 200[MeV], r//(n so) = 
C/( ra o-so) = 0.3 and KTo/(noSo) = 0.6. Note that uqSo represents the entropy density 
in the equilibrium state because sq is the entropy per particle number. 

As is expected, Fig.l shows that the Brillouin peaks owing to the sound mode 
is enhanced by the relativistic effects, while the Rayleigh peak owing to the thermal 
mode is the same as in the non-relativistic case. 
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2.2. In the case of particle frame 

Next let us take fluid dynamic equations in the particle frame. Does any differ- 
ence arise in the density spectral function depending on the choice of the frame. A 
typical equation in the particle frame is due to Eckart. 31 ) It is, however, well known 
that the Eckart equation shows a pathological behavior; 32 ) i.e, when the fluctuations 
around the thermal equilibrium is described by this equation, the fluctuation tends 
to diverge as t goes infinity; see Appendix B for a detailed discussion. 

Then, does any first-order relativistic fluid dynamic equation for a viscous fluid 
in the particle frame show such a pathological behavior? It is then noteworthy that 
a new fluid dynamic equation in the particle frame constructed from the relativistic 
Boltzmann equation by Tsumura, Kunihiro and Ohnishi (TKO) 5 ) turns out to be a 
stable one in the sense that any fluctuations around the thermal equilibrium state 
relax down to recover the equilibrium. 35 ) So let us take TKO equation in the particle 
frame and examine whether the spectral function of the density fluctuations shows 
any frame dependence. 

In the case of TKO equation in the particle frame, the dissipative terms are 
given by 

T » v = v [d^u v + - ^A^(d ± -u) } - C'(3u"u u - A^)(d ± -u) 

+ K{u^dlT + u u d^), (2-46) 
z/ 1 = , (2-47) 

where Q = C/(37 — 4) 2 with 7 being the ratio of the specific heats as before. 

From the same procedure as taken for the Landau equation, we obtain the 
spectral function as follows (the detailed derivation is given in Appendix C); 



with 



where 



{(6n(k,t = 0)) 2 ) v 7 ; w 2 + x 2 A; 4 

i [ (u-c s k) 2 + r£k* (cu + c s k) 2 + r£k* b { } 

Fb = \[X{1 - 1) + ^ TK ° " ^(«To + 3 C ')] ee r B TKO , (2-49) 
z riQCp 



V ™> ee (C + ±t,)/wo. (2-50) 

Note that since the fluid dynamic fluctuations around the equilibrium state is re- 
laxing in TKO equation (see Appendix C), we have obtained the spectral function 
without difficulty as was in the case of Landau equation. Moreover, the deviation of 
the width of the Brillouin peaks reads 

^ B TKO ^-^i( K T + 3C'), (2-51) 

ZUQCp 
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Fig. 2. The spectral function in full relativistic TKO equation (solid line) and its minimal relativis- 
tic version (dashed line) with the same parameters as those in Fig.l. Relativistic effects does 
not appear in the Rayleigh peak as in Landau case, but more enhance the Brillouin peaks than 
in the Landau case. Note that the scale of the vertical line is much bigger than that in Fig. 1. 



which is definitely negative. So the relativistic effect acts to enhance and sharpen 
the spectral function of the density fluctuation in the TKO equation in the particle 
frame in comparison with the minimal relativistic case than in the energy frame; see, 
however, the next subsection for the Israel-Stewart equation in the particle frame. 

If we estimate the relativistic effects using the EoS of massless classical ideal 
gas, we have 

^b TK ° = -^(^ + 3C7To). (2-52) 

It is noted that the effective bulk viscosity £ = C/(^7 — 4) 2 is finite in the massless 
case although 37 — 4 — > 0. 5 )' 35 ** The minimal relativistic part of the width is Eq.(2-44) 
with ( replaced by We see that the relativistic effect due to the thermal con- 
ductivity is the same as Landau equation but another relativistic effect due to the 
bulk viscosity exists. In addition the relativistic effect cancels out with some part 
of the minimal relativistic contribution due to the bulk viscosity and the thermal 
conductivity in Eq.(2-44) with £ replaced by (' . That is, the bulk viscosity and the 
thermal conductivity does not have any contribution in net in TKO equation for the 
massless classical ideal gas. 

The resultant spectral function is shown in Fig. 2 together with that of the min- 
imal relativistic case; the parameters are the same as in Fig.l: As in the case of the 
Landau equation, any relativistic effects does not appear in the Rayleigh peak while 
the Brillioun peaks are enhanced by them. The enhancement is more prominent than 
in the Landau case because the bulk viscosity is involved in addition to the thermal 
conductivity as relativistic effects. However, it is noteworthy that relativistic effects 
only appears in the Brillouin peaks but not in the Rayleigh peak, irrespective of the 
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choice of the frame of the fluid dynamic equation. 

2.3. In the case of Israel- Stewart equation in particle frame 

The Israel-Stewart (I-S) equation 33 ) is a second-order equation where relaxation 
times are contained and has a form of the telegrapher's equation. Thus in the 
I-S equation, the so called causality problem is formally resolved. However, the 
fluid dynamic equations are supposed to describes the behaviors of the conserved 
quantities in even larger scale than the mean free path, and the problem of causality 
takes place in short wavelength region. Therefore, it is expected that the acausal and 
causal fluid dynamic equations should give the same description for the phenomena 
with long wavelengths, or hydrodynamic modes. In fact we shall show that the I-S 
equation gives completely the same results for the spectral function of the density 
fluctuations as the Landau equation gives. 

In the I-S equation in the particle frame, the dissipative terms are given by 



t» v = -II A** + q p u u + q v u p + vr^, (2-53) 

and v p = 0. Here 

n = -C( V + Pou a d a n - a d^), (2-54) 

<f = KTA^(^d u T - u a d a u v - l3 lU a d a q v - a d v n + ai d a ir a u ), (2-55) 

tt^ = 2nA^ ab (d a u b - f3 2 u c d c 7r ab - ai d a q b ), (2-56) 



with u^q^ = 0, ir pv = ir up , u p tt^ v = and 7fy = 0. Here 0o, Pi and /% are the re- 
laxation time of the bulk viscosity, the heat flux and the shear viscosity, respectively, 
ao (cci) is the coupling of the bulk viscosity and the heat flux (the shear viscosity 
and the heat flux). A^ upa is a projector defined by 

A^ pa = -\A w A va + A^° ' A vp - -A^A pa ], (2-57) 
2 3 

Applying the similar procedure as in the first-order equations, we find the time- 
dependent density fluctuation is given by 

5n(k,t)/6n(k,0) ~ (1 - -)e~ xk ' 2t + - cos(c s kt)e- rBk2t 

7 7 

+0(k 2 ) x [e-W + e-'/ 2/32r > + e-wot/mwo-VW]] , (2-58) 

with 

= ~ 1)X + n + c 2 s T (k/w - 2 X a P )} = if, (2-59) 

which is exactly the same as given in the Landau equation in the energy frame. 
See Appendix D for a detailed derivation of these results. We note that the last term 
in Eq.(2-58) 

wot 



14 



Y.Minami and T.Kunihiro 



has an exponent which can be positive depending on the value of the relaxation time 
In fact this term is a remnant existing for the Eckart equation, which as shown 
in Appendix B, behaves pathologically as 

«5n(M)~exp[^-t]. (2-61) 

From Eq.(2-60), one can see that if the relaxation time of the heat current is large 
enough to satisfy the inequality 

Pi > — , (2-62) 
w 

the density fluctuation would relax to the equilibrium state. On the other hand, if 
the relaxation time of the heat current is so small as to satisfy the condition 

Pi < — , (2-63) 

the density fluctuation will not relax down and the system stays pathological even 
though the relaxation time is finite; that is, the Israel-Stewart equation in particle 
frame takes over the pathological behavior of Eckart equation. 

Here let us assume that the relaxation time is sufficiently large to satisfy the 
inequality Pi > Then the spectral function can be obtained as 

s nn (k,co) i 2 X k 2 i r r B k 2 



((Sn(k,t = 0)) 2 ) v 7 ; cj 2 + x 2 fc 4 7 L (w - c s k) 2 + r 2 k* 

+ i + o(**- 2/ « 



(u + c s ky + r 2 k 4> v ; l co 2 + i/(/3 C) 2 

yjhrj , 2w /[(/3iw - l)/cT ] 

+ W 2 + l/(2/3 2 7 ? )2 + ^ +u; 2 / p iU;o _ 1)Kro] 2J- ^ > 

Apparently, the spectral function has six peaks including the conventional three 
peaks, but the new three Lorentzian functions should vanish in the long wavelength 
limit k — > 0, because the strength of these is in the second order of k. Therefore 
Israel-Stewart equation gives completely the same result for the spectral function of 
the dynamical fluctuations as Landau equation does in the long wavelength limit, as 
shown in Fig.l. That is, the relaxation times does not affect the result in the fluid 
dynamical regime. 

Now one sees that the two relativistic fluid dynamic equations in the particle 
frame give different results on the dynamical density fluctuation; I-S equation gives 
the same result as that given by the equation in the energy frame whereas the TKO 
equation shows a frame dependence in the spectral peak of the sound mode. What 
is the origin of the difference? It may be attributed to that of the condition for 
the dissipative part of the energy-momentum tensor ST^ in these two equations 
though in the same particle frame: The I-S equation 33 ) is constructed with an ad- 
hoc postulate on 5T^ U as in the Eckart equation 31 ) that u^bT^Uy = which the 
Landau equation in the energy frame does satisfy because u^bT^ = by the very 
definition of the energy frame by Landau and Lifshitz. 30 ) On the other hand, the 
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TKO equation 5 ) in the particle frame is derived by a powerful reduction theory called 
the renormalization- group (RG) method 36 ) from the relativistic Boltzmann equation 
and it is found 5 )' 35 ) that their equation satisfies 5Tff = but not u^ST^Uu = 0. It 
is to be noted that the RG method gives the Landau equation for the fluid dynamic 
equation in the energy frame. 5 ) It is clear that more work is needed to establish the 
correct fluid dynamic equation in the particle frame. 



§3. The behavior around the QCD critical point 



We analyze the behavior of the spectral function of the density fluctuations 
around the QCD CP on the basis of the dynamic as well as static scaling laws. 22 )' 24 ) 
We introduce the static critical exponents 7 and a which are defined as follows 

c n = c t~ & , K T = K t-\ (3-1) 

where t = \(T — T c )/T c \ is a reduced temperature, cq and Kq are constants and 
Kt = (l/no)(dn/dP)T is the isothermal compressibility. Using the known thermo- 
dynamical identities 22 ) 

dP\ 2 



c P = c n + (To/no) l^—j K T , (3-2) 

cl = -4^—, (3-3) 
w c n K T 

ap = KT ' (34) 

we see that the specific heat at constant pressure, the sound velocity and the thermal 
expansivity at constant pressure has the critical behavior as 

cp ~-^UfJ n ' ' (3 - 5) 

«p~^o(J0 t~\ (3-7) 

respectively. Since the phase transition at the QCD CP belongs to the same univer- 
sality class Z2 as the liquid-gas CP, the critical exponents are given as a ~ 0.11 and 
7 ~ 1.2. Recently, it has been argued 20 ) on the basis of the analysis of the liquid-gas 
CP by Onuki 37 ) that the bulk viscosity may show a singular behavior around the 
QCD CP. So we introduce the critical exponent as 

C = Coi" ac , (3-8) 

where Co is constant. The exponent for the liquid-gas CP is predicted to be 
zv — a. 37 ) Here z is the dynamical critical exponent and v is the exponent which 
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represents the singularity of the correlation length. These are given by z ~ 3 and 
v ~ 0.63 around the liquid-gas CP. 

We denote the exponent of the thermal conductivity by a K , 

K = K t~ aK , (3-9) 

where kq is a constant. It is known that a K ~ 0.63 around the liquid-gas CP. Let us 
take it for granted that the exponents and a K at the QCD CP are given by those 
at the liquid- gas CP. That is, we take 

zv — a ~ 1.8, (3-10) 
a K ~0.63. (3-11) 

Combining these ingredients including the possible singular behavior of (, we see 
that the width of the Rayleigh peak behaves as 

r R ~^(|^)Vv (3 . 12) 

K Q T \dPj n { 1 

Then, we see 

r R ~^-°-, (3-13) 

which shows that the width Fr becomes narrow as the QCD CP is approached. 
We emphasize that this result is independent of the choice of the relativistic fluid 
dynamic equation or frame. 

So far for the width of the Rayleigh peak. How about the width Fq of the 
Brillouin peaks? For the Landau and I-S cases, it has the critical behavior as follows, 

r B ~— (3-14) 

We note that this singularity comes from that of the bulk viscosity. 

In the case of TKO equation, we first note that the critical behavior of the 
effective bulk viscosity is given by 

which shows that the effective bulk viscosity has a positive exponent and does not 
show a singular behavior because ~ 1.8 and 7 ~ 1.2. Instead, the singularity of 
the Brillouin peaks for TKO equation comes from that of the thermal conductivity; 

r B ~r^ KO t-^- & \ (3-16) 

with 

™o \dTJ n _ 



j-iTKO _ K 



2c n 

We note that the strength of the divergence of Fq for TKO equation is weaker than 
that for the Landau and I-S equations. 
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Fig. 3. The spectral function at t = 0.5 and k = 0.1 [1/fm]. The solid line represents the 
Landau/Israel-Stewart case, while the dashed line the TKO case. The strength of the Bril- 
louin peaks becomes small due to the singularity of the ratio of specific heats. 
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Fig. 4. The spectral function at t = 0.1 and k = 0.1 [1/fm]. We see that the Brillouin peaks 
which correspond to sound wave dies out and the difference between the Landau and TKO 
cases disappears. Note that the scale of the vertical line is much bigger than that of Fig. 3. 



Anyway, we have confirmed that the width Fq may diverge at the QCD CP 
irrespective of the relativistic fluid dynamic equations. 

Unfortunately or fortunately, these singular behaviors of the width of the Bril- 
louin peaks around the QCD CP may not be observed: Note that the strengths of 
the Rayleigh and the Brillouin peaks are given in terms of 7, the ratio of the specific 
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heats, which behaves like 

7 = Cp/cn ~ t~ x < +& -> oo, (3-18) 

in the critical region. Then the strength of the Brillouin peaks is attenuated and 
only the Rayleigh peak stands out in the critical region, as follows; 

S nn (k,u) ~ ((5n(k,t = 0)?) J^l k , , (t~0). (3-19) 

Figure 3 and 4 show how the spectral function behaves around the QCD CP for 
k = 0.1 [1/fm], t = 0.5 and t = 0.1, respectively. We see that the strength of the 
Brillouin peaks becomes small and dies out as the system approach the QCD CP. In 
addition, the static correlation function ((6n(k,t = 0)) 2 ) has the critical behavior as 

((5n(k,t = 0)f)~ T ^L ¥ , (3-20) 

with B = TKt/V, where V is the volume. 38 ) Equation (3-20) would show a singular 
peak in the forward direction k = 0, because Kt ~ i" 1 ' 2 and £ 2 ~ i~ L26 . *) Then 
the strength of the Rayleigh peak will be most drastically enhanced in the forward 
angle in the critical region. 

The critical behavior on the spectral function of the density fluctuation is sum- 
marized as follows. 

1. Disappearance of the strength of the sound mode: 

The sound mode loses its strength at the QCD CP because of the divergence 
of the ratio of the specific heats. Thus, we see that the formation of the sound 
mode itself is suppressed. 

2. Divergence of the width of the sound mode Fr: 

According to Eq.(3T4), the width of the sound mode diverges at the QCD CP, 
which means that the sound mode is strongly damped and attenuated in the 
critical region. Therefore we see that the sound mode is very much attenuated 
even if it is created with a small strength. 

3. Narrowing of the width of the thermal mode Tr: 

The width of the thermal mode vanishes at the QCD CP because of the diver- 
gence of the specific heat at constant pressure. This means that the entropy 
density exhibits critical slowing down and the thermal mode is the critical mode 
at the QCD CP. That is, the entropy density is the order parameter and mainly 
governs the critical dynamics at the QCD CP as in the liquid-gas CP. 

§4. Discussions 

We have seen that the sound mode would lose its strength and be strongly 
damped as the system approaches the critical point. But why at all do sound modes 
die out at the critical point? To answer this question, let £ = iot~ v be the correlation 

This is known as the critical opalescence. 22 '' 26 ' 
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length which diverges as the critical point is approached. If we write the wave length 
of the sound mode by A s , the fluid dynamic regime is expressed as 

i « A s , (4-1) 

with which condition the sound mode can develop. 29 ) However, in the vicinity of 
the critical point, the correlation length £ becomes very large and eventually goes to 
infinity, so the above inequality can not be satisfied, and the sound mode can not be 
developed in the vicinity of the critical point. 29 ) 

From this argument, we can speculate that phenomena inherently related to the 
existence of the sound mode may disappear around the critical point. One of such 
phenomena is the possible Mach cone formation 39 ) by the particle passing through 
the medium with a speed larger than the sound velocity c s . Such a Mach-cone 
like particle correlations are observed in the RHIC experiment. 40 ) If such three- 
particle correlations have been confirmed to be a Mach-cone formation, then the 
disappearance or suppression of the Mach cone would be a signal that the created 
matter has passed through the critical region, showing the existence of the QCD 
critical point. Even if the thermal wake also contributes to the formation of Mach 
cone, 41 ) a suppression of Mach cone may be expected by the attenuation of the sound 
mode. So it would be very interesting to see possible variation of the strength of 
the Mach cone according to the variation of the incident energy of the heavy-ion 
collisions. In theoretical side, it is an intriguing task to explore the fate of Mach 
cone with an EoS which admits the existence of the CP. 

Precisely speaking, the analysis based on fluid dynamics is, however, only valid in 
the fluid dynamical regime: k£ « 1, that is, our discussions on the critical behavior 
are based on extrapolation from the fluid dynamical to the critical region: k£ » 
l. 29 ) This extrapolation is known to make good sense for the non-relativistic case. 22 ) 
It is well known that the coupling of the thermal mode to the diffusive transverse 
mode becomes essential for the description of the dynamical critical phenomena in 
the close vicinity of the CP for the non-relativistic case. 23 ) This coupling can be 
analyzed by the mode-mode coupling theory. 23 )' 24 ) Such an analysis is interesting 
but beyond the scope of the present work. 

§5. Summary and concluding remarks 

We have explored how the singularities of the thermodynamic values as well as 
the transport coefficients affect the dynamical density fluctuations around the QCD 
critical point(CP) using dissipative relativistic fluid dynamics. Our analysis is an 
extension of those made for non-relativistic case to the relativistic case. We have 
shown that the sound modes which are directly related to the density fluctuation are 
enhanced by relativistic effects, but tend to be attenuated around and would eventu- 
ally die out at the CP and, hence, are not the soft mode. Our analysis based on the 
relativistic fluid dynamics has shown that the genuine and remaining soft mode at 
the QCD CP is the diffusive thermal mode. Our analysis also suggests that the possi- 
ble divergent behavior of the bulk viscosity may not be observed through the density 
fluctuations. We have also argued that the possible suppression or disappearance of 
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Mach cone can be used as a signal of the existence of the QCD CP. 

We have shown that the relativistic effect on the spectral function of the density 
fluctuation appears only in the width of the Lorentzian function. More precisely, 
the relativistic effect enhances only the Brillouin peaks due to the sound mode but 
not the Rayleigh peak irrespectively of the equations and the frames. Therefore one 
might tend to think that the choice of the frame is just a matter of convenience and 
gives no physical difference. However, it is noteworthy that there is a slight difference 
in the enhancement depending on the choice of the fluid dynamic equation in the 
particle frame. The spectral function obtained in the Israel-Stewart(I-S) equation 
gives the same result as that in the Landau equation in the energy frame whereas 
the equation by Tsumura-Kunihiro-Ohnishi(TKO) 5 ) gives a more enhanced peaks to 
the sound mode than in the two equations. We have suggested that the origin of the 
difference may be attributed to that of the condition for the dissipative part of the 
energy-momentum tensor 8T^ V in these two equations though in the same particle 
frame. It is certain that more work is needed to establish the correct fluid dynamic 
equation in the particle frame. 42 ) 

We have also elucidated the properties of some of dissipative relativistic fluid 
dynamic equations: In the case of Eckart equation, the density fluctuation around 
the equilibrium is unstable and tends to diverge as t goes infinity. This pathological 
behavior is due to the heat current induced by the time derivative of the fluid veloc- 
ity. Since this property is carried over to the Israel-Stewart(I-S) equation in particle 
frame, the density fluctuation in the I-S equation also behaves pathologically if the 
relaxation time of the system happens to be small, although it is formally causal; 
if the relaxation time is large enough, the I-S equation in particle frame gives com- 
pletely the same spectral function for the density fluctuations as Landau equation 
gives. That is, formally acausal fluid dynamic equations suffices in description of 
phenomena in the fluid dynamical regime with long wave lengths if the equation 
gives a relaxation of the fluctuations around the thermal equilibrium. 

In the present work, we have confined ourselves to the study of fluctuations 
around the thermal equilibrium. To make a direct relevance to the relativistic 
heavy-ion collisions, study on the fluctuations around expanding flow background 
like Bjorken flow should be done, which constitute one of the future works. It should 
be also an intriguing task to explore the fate of Mach cones or shock waves in general 
in the critical region by explicit calculations. 

It is important to study the coupling between the thermal fluctuations and the 
transverse mode using the mode-mode coupling theory in the close vicinity of the 
critical region for the relativistic case. The present work constitutes the basis for 
such a more complete analysis. 
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Appendix A 

derivation of thermodynamic identities 

In this Appendix, we derive thermodynamic identities we have used. Here s is 
not the entropy density but the entropy per the particle number, c n = Tq(8s / dT) n 
and cp = To(ds/dT)p are the specific heats at constant density and pressure, re- 

1/2 

spectively, c s = (dP/de) s is the sound velocity, ap = —(l/no)(dn/dT)p is the 
thermal expansivity at constant pressure, and 7 = c~p/c n is the ratio of the specific 
heats. §|^y represents Jacobian determinant: 

d(a, b) _ da db da db 

9M = dc~dd ~ dd'dl' ( ' ' 

A.l. Derivation of (dP/dn)T = {vjQC 2 s )/{n^) 

(-) 



where 



d(P,T) 
d(n,T)' 
d{P,T) d{P,s) d(e,s) d{n,s) 
d(P,s) d(e,s) d{n,s)d{n,TY 




(A-2) 



d(e,s) 
d(n, s) ' 
d(e, s) d(e, n) 
d(e, n) d(n, s) ' 




(A-3) 
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Using the first law of thermodynamics 

de = Td(ns) + fj,dn, 



we find 



and 



96 ^ =T^( d 4\ =nT, 



d(ns) J n \ds 
d(ns)\ fi 



s + n 



9n y e T' 
f ds\ \i 



\dnj e T' 



where we have used e = Tns — P + fin and u; = e + P. 
Then we obtain 

f de\ w 
\dn) s n 

and so 

' dP\ w c 2 



dn J T no7 
A. 2. Derivation of (dP/dT) n = woapc^/j 

dP\ _d(P,n) 



drj n d(T, n y 

__ d(P,n) d(T,P) d(s,P) d(s,e) d(s,n) 
~ d(T,P) d{s,P) d{s,e) d(s,n) d(T,n)' 
dn\ ( dT\ ( dP\ ( de\ ( ds 



dT J p \ds J P \de J s \dnJ s \dT 



T 2 w 
noctp ■ — • c • — 
c P n 

WQOipC 2 



7 

A. 3. Derivation of(ds/dn)T = — (wQCtpc 2 s ) / {n^) 

Assuming that the total particle number N is constant, we find 

ds_\ ( d{S/N) \ 
dnj T \d{N/V)) T ' 

1 dV (dS\ 
~ Wd{l/V) \dV) T ' 

n 2 Q \dV. 
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Using Maxwell relations 

dS\ (dP\ 



ov) t -\ot) v < a -»> 



We obtain 

ds \ 1 (dP\ 



dn J rp n"l\dT ) y 
1 fdP s 

ds\ WQCtpc 2 s 



dn J rp 7Jq7 
where we have used Eq.(A-8). 

A. 4. Derivation of (woToapc 2 ) / (nocpj) = 1 — 1/7 



Cp — lQ 



= T 



d(T, P) ' 
8(s,P) d(T,n) 
d(T, n) d(T, P) ' 

/ ds\ f dn 



c P 



ds 


) +T 

/ n 


df 


ds 


\ f dn 


dn 


JA&T 




WQCtpC 2 s 


(" 





\dn J T \ dT J p 



Dividing the both side by cp, we obtain 



(A-12) 



p 

= T -[-^P)-(-n ap). (A-13) 



l-i = ^M. (A-14) 

7 n cp7 

Appendix B 

In the case of E chart equation (-particle frame) 



In this Appendix, we shall show explicitly how the pathological behavior of the 
Eckart equation hinders the procedure of the calculation of the spectral function of 
the density fluctuations for an instructive purpose. 

In Eckart equation, the dissipative terms are given by 



T ^ = rj [d^u u + d u ± u^ - -A^(d ± -u)} + (A» u {d ± -u), 



2 

+ K^d^T + u v &\T - 7V> • d)u v - Tu v {u ■ (B-l) 
= . (B-2) 
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The term —Tu^(u ■ d)u v — Tu u (u ■ d)u^ in Eq.(BT) represents the heat current 
induced by the time derivative of the fluid velocity, and this term is identified to the 
origin of the pathological behavior. 

The linearized fluid dynamic equation reads 



d5n 



d5v 



+ n V • Sv= 0, 

d5T y 



w - 



dt 



v V 2 5v - (-r, + C )V(V • 5v) + V(6P) - «V(— ) 

0, 



d 2 5v 



no- 



85s 
~dt 



reV 



85v 



—V 2 5T= 0. 
To 



(B-3) 
(B-4) 

(B-5) 
(B-6) 



It is noted that the second order derivative in time appear in Eq.(B-5) which arises 
from the heat current induced by the time derivative of the fluid velocity. From the 
same procedure as was done for the Landau equation in § 2, we obtain 



(5n(k, z) 
Sv\\(k,z) 
5T{k,z), 



Sn(k,0) \ 
0) + (-z^ + l)S V]l - ik^8T(k, 0) 



"wq dt 
ape 2 



w 



-^Sn(k, 0) - ik^-fy (k, 0) + ^5T(k, 0) J 



(B-7) 



where 



A 



ikno 



ih^f- -z 2 ^ + » t k 2 

u>o7 wo 1 

^ -izk-*- 

no 7 wo 



— Z 








The Fourier-Laplace transform of the density fluctuation is given by 



(B-8) 



5n(k, z) = {(A _1 )n - — ^-(A _1 )i 3 }(5n(fc, 0) 



a P c 2 s 



n 7 



(B-9) 



where the irrelevant terms which do not contribute to the density fluctuations are 
omitted, det A is evaluated to be 



det A 



c„n 0r a kT 3 2 kT 
1-2* 1- - Ar(7x h 



o 



• - , • - )}: 

T w w w n woc n 

+z 2 k 2 ( 1X + 2xa P c 2 s T ) + z(k 2 c 2 + k^xn) + k\c 2 s }. 



(B-10) 



The fourth order term in z appears in det A because of the heat current induced by 
the time derivative of the fluid velocity. We can factorize this to second order in k 



KC n no w 2 

del .1 ~ 



Wn 



(z - -^){z + X k z ){z + r B k 2 + ic s k)(z + r B k z - ic s k), (B-ll) 

K±o 



where 



= ^[(7 " 1)X + n + ^T (k> - 2 X a P )] = r| c . 



(B-12) 
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Comparing Eq.(BTl) with the case of Landau equation, we see a new mode is 
present in the present case, which is found to be an unstable mode so that the 
density fluctuation increases with t without limit, 



Sn(k,t) ~ exp[^-t]. 



(B-13) 



Of course, this is an unwanted result because any fluctuation around the equilibrium 
state should relax down to recover the equilibrium state. This pathological properties 
of the Eckart equation was first noted by Hiscock and Lindblom. 32 ) 



Appendix C 

Detailed derivation in TKO equation in particle frame 



The linearized TKO equation 5 ) reads 



d5n 
Ht 



+ uqV ■ Sv = 0, 



w o^r ~ ?/V 2 (5i> - (-rj + C')V(V • Sv) + V(5P) - «V(^^) = 0, 



dt 



dSs 3C 



dt 

dSv K n rm 



Accordingly, we have 

(Sn(k, z) 
5v\\(k,z) 
ST(k,z) t 



a P c% 



Sn(k,0) 
6v\\(k,0)-ik^6T(k,0) 



WO 

^5n(fc,0)-ifc 5 g s fo||(fc,0) + 7W 



where 



A = 



ikno 



\ 



ik-^ 3 - 

"07 
a P c 2 „ 



zkK 



(C-l) 
(C-2) 

(C-3) 



(C-4) 



(C-5) 



Using these result, we have the density fluctuation up to the second order in k, as 
follows, 

„2 



5n(k,z)/5n(k,0) = [(A _1 ) n 



n 7 



(A 



)13\ 



(z + r B k 2 + ic s k)(z + r B k 2 - ic s k) - ±cjk 2 + z\k 
(z + xk 2 )(z + r B k 2 + Ac s k){z + r B k 2 - ic s k) 



Here 



with 



r B = - 2 lxh - i) + ^ TKO 



ape 2 
n c P 



(kT + 3( )]=T£ 



TKO 



v? K ° = (C + iv)/w , 



- (C-6) 

(C-7) 
(C-8) 



26 



Y.Minami and T.Kunihiro 



Thus the spectral function is now given by 



Snn{k, < 



((<5n(M = 0)) 2 ) 



(1 



1 



2 X k 2 



1 

+ - 



7 cj 2 + x 2 k 4 ' 

r B k 2 



r B k 2 



j\cj-c s k) 2 + T 2 k 4 (lo + c s k) 2 + r 2 k 4 > 

Appendix D 

Detailed derivation in the case of I-S equation 



(C-9) 



Applying the similar procedure as in the first-order equations, the linearized I-S 
equation reads; 



d5n 
~dt 



+ n-oV ■ 8v = 0, 



(eo + Po) ^ + V(<5P + ,577) + ^ + V ■ for = 0, 

or 7-7- 

<5i7 + C[V • foj + Po— - a V • 5q] = 0, 

<5g - KT [-^rV(5T) - ^ - - a V(<5iT) - ftl V • for] = 0, 

• • • 2 • • i9for iJ 

(5tt^ - ril&Sv 3 + dW - -g^V ■ 5v - 

-a 1 (d i 5q j + &> 5q i - -g^V • Sq)] = 0. 
o 



(D-l) 
(D-2) 
(D-3) 
(D-4) 
(D-5) 

(D-6) 



Performing Fourier-Laplace transformation, we find for the longitudinal compo- 
nents, 



A 



( 5n(k,z) \ 
Sv»(k, z) 
5T(k,z) 

sn(k,z) 

Sq\\(k,z) 

\57Tu n(k, Z )J 



( 6n(k,0) \ 

(3oCSn(k,0) 
-(c 2 s ap/n "/)5n(k,0) + (n c n /woT )5T(k, 0) 
kToSv^O) +/3iKT ^||(fc,0) 
6v\\{k,0) + (l/w )6q\\(k,0) 

2^/32^1111(^,0) / 



(D-7) 



where dq\\ 




k ■ 5q, for || y = 


k l -k> ■ 


TTij and 








( 


z 


ikno 











\ 









ik( 





1 + zM 


— ikao( 





A = 




no 7 







-v n c n 
w T 

kn 




kctQKTo 


ik 
wo To 

1 + z/?ikT 




ka\ 




V 


k^- 

no7 




z 


k 7 



j_ 

wo 




z 

wo 

— pair] 


_k_ 

wo 

1 + 2/W 



(D-8) 
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For simplicity, we have neglected the coupling terms ao and a±, which are expected 
small and irrelevant for the following argument. 

Applying the similar argument done in the case of Landau equation, we have 



5n(k, z) = [{A~ l )u - C -^L{A- 1 ) n ]5n{k, 0). (D-9) 

no7 

In the present case, det A is given by 

det A = F lZ 6 + F 2 z 5 + F 3 z A + F±z 3 +k 2 (F 5 z A + F 6 z 3 + F 7 z 2 + F 8 z) 

+k\F 9 z 2 + F 10 z + F n )], (D-10) 

where F\ ~ Fu are defined by 

Fi = 2f3 (3 2 ( V KT ((3 lW o - 1), (D-ll) 

F 2 = 2 /W? + «T (Awo - 1)(A>C + 2/321?), (D-12) 

F 3 = kT ((] iW0 - 1) + ^o(/?oC + 202V), (D-13) 

F 4 = u> , (D-14) 

F 5 = /3iCr/AcT (^/3o + 2/3 2 ) + 2/WW7X - 2xc 2 apT + ^kT q ), (D-15) 

F 6 = PxkTqWqVi + io (A) 

+2/3 2 r/)( 7 x - 2 XC 2apTo + c^ lK T ) 

+Cr ? (^/3o + 2/? 2 + 2 WoC 2/3 /3 2 ), (D-16) 

^7 = ^o[7X - 2x^ap3b + c^Tq + ^ + c 2 (/? + 2/3 2V )}, (D-17) 

= ™oc 2 , (D-18) 

^9 = 7X(A)C + 2/32»/) + 2(3 p 2 ( V w c 2 sX , (D-19) 

*io = «>obw + c sX(A)C + 2/3 2 r7)], (D-20) 

*n = ^X, (D-21) 



which are all independently of A; and z, As before, det A can be approximately 
f actor ized for small k, 

det A ~ ^(/3oC^ + 1 + O(fc 2 ))[(/3iu;o - 1)«T + u; + 0(fc 2 )] 
x(2/3 2 r ? z + l + 0(/c 2 ))(z + xfc 2 ) 

x (z + r B /c 2 + ic s k)(z + T B A; 2 - ic s k), (D-22) 

with 

T B = ^[(7 - 1)X + ^ + c 2 T (k>o - 2 X a P )] = /f, (D-23) 

which will be found to be the width of the Brillouin peaks and is equal to that for 
the Eckart equation, r® c . One can see that the relativistic effect enter two ways; one 
is due to thermal conductivity which tends to enhance the width while the thermal 
expansivity would reduce the width. The fact that the relativistic effects might 
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possibly enhance the Brillouin width is in contrast to the case of TKO equation 
where the relativistic effects only tends to reduce the width of the sound modes. As 
we will see in the text, the net relativistic contribution to .T B S is found to be negative. 
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